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ABSTRACT 

We present a high quahty CMD for a 36' x 36' field located 8° (~ 7 kpc) from 
the LMC center, as well as a precise determination of the LMC surface brightness 
derived from the resolved stellar population out to this large galactocentric radius. 
This deep CMD shows for the first time the detailed age distribution at this 
position, where the surface brightness is \^ ^ 26.5 mag/D". At a radius 474' 
the main sequence is well populated from the oldest turnoff at / ~ 21.5 to the 
2.5 Gyr turnoff at / ~ 19.5. Beyond this radius, a relatively strong gradient 
in the density of stars with ages in the ~ 2.5-4 Gyr range is apparent. There 
are some stars brighter and bluer than the main population, quite uniformly 
distributed over the whole area surveyed, which are well matched by a 1.5 Gyr 
isochrone and may be indicative of a relatively recent star formation, or merger, 
event. The surface brightness profile of the LMC remains exponential to this 
large galactocentric radius and shows no evidence of disk truncation. Combining 
the information on surface brightness and stellar population we conclude that the 
LMC disk extends (and dominates over a possible stellar halo) out to a distance 
of at least 7 kpc. These results confirm that the absence of blue stars in the 
relatively shallow off-center CMDs of dirr galaxies is not necessarily evidence for 
an exclusively old stellar population resembfing the halo of the Milky Way. 

Subject headings: galaxies: Magellanic Clouds — galaxies: halos — galaxies: evolution — 
galaxies: stellar content 

1. Introduction 

The nature of the stellar populations at the outer edges of dwarf irregular (dIrr) galaxies, 
and in particular whether these galaxies have a truly old stellar halo has been a subject of 
some controversy over the last few years. The presence of a background sheet of red stars, 
extending beyond the most prominent blue, young stellar population in 1C1613 and other 
dIrr galaxies was discussed already by Baade (Sandage 1971). But, what is the nature of this 
sheet of red stars? Minniti & Zijlstra (1996) claimed that the Local Group dlrr WLM has 
an old (as old as the Milky Way globular clusters) stellar halo, from the sole evidence that 
the color-magnitude diagram (CMD) of an outer field extending out to 8' from the center of 
the galaxy, along its minor axis, does not contain blue stars down to I ~ 23 (M/ ~ —1-8), 
and shows a well defined tip of the RGB. However, since all stars older than about 1 Gyr 
produce a RGB at very approximately the same location, and main-sequence (MS) stars 
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older than a few hundred Myr would be below the hmit of Minniti & Zijlstra's photometry, 
all that can be confidently concluded from the lack of bright MS stars is that the stellar 
population is older than a few hundred Myr (see also the discussion by Aparicio, Tikhonov 
& Karachentsev 2000). Other evidence, such as the presence of C stars in the outer regions 
of dwarf galaxies (e.g. Letarte et al. 2002), or a full modeling of the CMD using synthetic 
CMDs (Aparicio & Tikhonov 2000) suggests in fact the presence of a range of ages. 

Determining the age structure and kinematic properties of the outer red component 
of dlrr galaxies may hold important clues on the process of dwarf galaxy formation and 
evolution. Do these stars represent a stellar component truly differentiated from the younger 
one, with the characteristics of a halo population similar to that of the Milky Way (old, metal 
poor, devoid of gas, kinematically hot), or is there a smooth gradient of stellar ages, which 
might suggest that star formation has been "shrinking" with time? (Hidalgo, Marm-Pranch 
& Aparicio 2003) 

A complete answer to these questions requires both very deep photometry, reaching 
the oldest MS turnoffs, and relatively high dispersion spectroscopy to obtain kinematics 
of individual stars. Prom deep CMDs at different galactocentric distances it is possible to 
characterize the ages of the stellar populations present, and to assess, therefore, whether a 
"break" in the stellar population, indicating an edge to the zone of protracted star formation 
and a transition to an old stellar halo actually occurs. The observational requirements for 
this task are, however, challenging. For the Milky Way satellites, it is necessary to reach 
1/ ~ 24 to obtain good precision photometry of the old MS turnoff region. For the rest of the 
Local Group, y ~ 28 is necessary. While the first is attainable from the ground, the latter 
is only possible using the HST equipped with ACS. Obtaining the required spectroscopic 
information is a highly challenging enterprise. 

The Magellanic Clouds are the two nearest Irr galaxies. Their old MS turnoffs are 
located near y ~ 23. While this magnitude is easily attainable from the ground, the huge 
angular sizes of the Clouds makes the mapping of their full extent and the characterization 
of their stellar population gradients a daunting task. As discussed by Irwin (1991) on the 
basis of isopleth maps, the LMC extends over an eUiptical area covering 23° by 17°, with the 
major axis oriented roughly in the North-South direction. This is also basically the extent 
of the LMC cluster system. On the North side. North of declination -62° (~ 8° from the 
center), the isopleths abruptly separate from each other. Does this indicate some kind of 
break in the characteristics of the stellar populations? Some information about the ages of 
the stars present at this position comes from the study of Kunkel, Irwin & Demers (1997), 
who searched for C stars in the periphery of the two Magellanic Clouds. Interestingly, they 
found a relatively sharp boundary for the outermost distribution of C stars at a radius of ~ 8° 
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along the major axis, which may indicate the actual maximum extent of an intermediate- age 
population. A detailed characterization of such a population, however, requires obtaining 
CMDs reaching the oldest MS turnoffs. A full dynamical study of the C star population of 
the LMC based on 1043 objects out to 8° from the center was published by van der Marel 
et al. (2002), who inferred the presence of a dark- matter halo containing roughly 2/3 of the 
total mass of the LMC. The presence of a kinematically hot, metal-poor, old stellar halo has 
been inferred by Minniti et al. (2003) from the kinematics of a sample RR Lyrae stars. 

With the purpose of characterizing the stellar population gradients in the LMC, we 
imaged four fields situated at distances ranging from 3.1 to 8 degrees from the LMC center 
using the CTIO 4m MOSAIC wide field imager. The full analysis of the whole dataset is 
underway (Gallart et al. 2004). In this Letter we analyze the CMD of the outermost field, 
which shows a somewhat surprising result on the age composition of the outer LMC, and 
use the whole dataset to trace the surface brightness across this galaxy. 

2. Observations and data reduction. 

We obtained V and I images of a LMC field centered at 7°8 from its center to the North 
{a : 05:13:17; S : -61:58:40) on January 17th and 18th, 2001, using the Mosaic II CCD Imager 
on the CTIO Blanco 4m telescope. The Mosaic II provides a total field of 36' x 36', with a 
resolution of 0'!27/pixel. Seeing was typically around l'(0. Exposure times were 6 x 900 s in 
V plus 8 X 600 s in I. Several short exposures were also obtained to retrieve the photometry of 
the brightest stars. The nights were photometric, and the Landolt (1992) standard star fields 
SA92, SA95, SA98, SA104 were observed several times each night for calibration purposes. 
We used Stetson's (2000) photometry for these fields to calibrate our photometry. During 
the same run and in another one in December 1999, we obtained photometry of similar 
quahty for three more internal fields. These results will be described elsewhere (Gallart et 
al. 2004), and used occasionally here (Section 3). 

The Mosaic frames were reduced in a standard way, as indicated in the CTIO Mosaic 
II Imager User Manual (Schommer et al. 2000), using the MSCRED package within IRAF. 
Finally, profile-fitting photometry of the LMC fields was obtained using the DAOPHOT / 
ALLFRAME suite of computer programs (Stetson 1987, 1994). A full description of the 
data reduction and photometry is provided in Gallart et al. (2004). 
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3. The LMC surface brightness profile: reaching the outer edge? 

We have derived the surface brightness for the four fields we observed in the LMC, 
by integrating the fiux in the CMD^. This method has the advantage that a very accurate 
field subtraction can be performed at the very low surface-brightness regime as foreground 
stars and background galaxies are for the most part well separated from the characteristic 
sequences of the LMC CMD. This is particularly important in the outermost field, where the 
flux contribution of Milky Way foreground objects exceeds by far that of the galaxy. The 
process will be fully described in Gallart et al. (2004). In short, we used the predictions of 
the Besangon Galaxy model (Robin et al. 2003), but scaling the number of predicted Milky 
Way stars to the number observed in our field in zones in the CMD clearly devoid of LMC 
stars. Without this scaling, we found that Besangon models predicted 25% more Galactic 
stars than we observed at relatively bright magnitudes. In addition, we only considered 
stars fainter than V — 15 when integrating the flux in the fleld, since no stars belonging 
to the LMC were observed brighter than this hmit. Since fleld stars brighter than V — 15 
are the main contributor to the total flux of the fleld, this greatly reduced the error in its 
subtraction. We neglected the contribution of background unresolved galaxies, since they 
are faint. 

In Figure 1, the surface brightness of our flelds measured in V (dots) is represented as 
a function of distance from the center of mass of the LMC. We have divided the outermost 
fleld into seven sections at increasing distances from the center of the LMC and calculated 
the surface brightness for each, in order to check whether the surface brightness variation 
remained smooth (which seems to be the case) or whether some kind of truncation was 
observed instead. The error bars were estimated by repeating the calculations with sky 
backgrounds that varied by 25% from our best estimate. We have added (triangles) two 
measurements from Hardy (1978). Exponential laws for the LMC disk as calculated from 
the whole dataset (Hardy plus our data, solid fine), and as derived by Hardy (long-dashed 
fine) and Bothun & Thompson (1988, short-dashed line) are represented. Our flt has a slope 
similar to that given by Bothun & Thompson, but is brighter. We also don't see evidence 
of the truncation found by these authors at a galactocentric distance of approximately 170'. 
We conclude that the surface-brightness measurements displayed are consistent with an 
exponential law out to the last point measured. 



^Thc contribution to the total flux of the unresolved LMC population (i.e. by the stars fainter than 
My — 5) may range between 1 and 10% of the total flux, as estimated from synthetic CMDs for a composite 
population (constant star formation rate from 15 Gyr ago to the present time) and for a purely old popula- 
tion. We estimate, therefore, that we may be underestimating the surface brightness by a maximum of 0.1 
magnitudes. 
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The field we will discuss in this Letter is the outermost one, situated from 450' to 490' 
from the center of the LMC. The surface brightness in this field ranges from 26.2 to 26.8 
mag/n"in V and is one magnitude brighter in /. To put the location of our field in context, 
we have drawn an arrow at the location of the outermost radius for which the LMC surface 
brightness was measured by Bothun & Thompson (1988). The tip of the arrow indicates an 
estimate of the V magnitude measured by these authors at this position, calculated using 
their B magnitude and assuming a color {B — V) ~ 0.5. This location also approximately 
coincides with the outer radius of the HI measured by Kim et al. (1998). The radius of 
the outermost field that we have measured is roughly equal to the radii of the maximum 
extension of the LMC cluster system (Olszewski et al. 1988) and the outer edge of the C-star 
system obtained by Kunkel et al. (1997). 



4. The Color-Magnitude Diagrams 

Figure 2 shows [{V — /),/] CMDs of the LMC, in order of increasing distance, R, from 
the LMC center, from galactocentric radius 450 to 490' (see caption for details). Stars with 
at least one valid measurement in each filter have been selected using the following limits 

for the error and shape parameters given by ALLFRAME: a(^v-i){= \/<Ty + aj) < 0.1, 
I sharp \< 0.3 and x ^ 4.5. A total of 55000 stars down to / < 24.0 have been measured. 
The CMDs reach a couple of magnitudes below the old MS turnofi' and therefore, this point 
is always cleanly reached with good photometric accuracy. Isochrones from Pietrinferni et 
al. (2004) for Z^O.OOl and ages 13.5 and 8 Gyr, Z=0.002 and age 4.5 Gyr and Z=0.004 and 
ages 2.5 and 1.5 Gyr have been superimposed. 

These CMDs show for the first time the details of the age structure of the stellar pop- 
ulation at the outer edge of the LMC. Out to a radius R~474' the MS appears populated 
quite smoothly from the oldest turnoff at / ~ 21.5 to the 2.5 Gyr turnoff at / ~ 19.5. 
Beyond this radius, a relatively strong gradient in the density of stars of ages ~ 4-2.5 Gyr 
begins to be apparent, with very few stars in this range beyond R=480'. There is a number 
of stars brighter and bluer than the main population, which are well matched by a 1.5 Gyr 
isochronc. They are distributed quite uniformly over the whole field. No foreground stars 
or background galaxies, however, are expected to populate the CMD at this position, and 
therefore they must be young stars genuinely belonging to the LMC. They may indicate a 
small burst of star formation 1.5 Gyr ago, the spatial extent of which cannot be evaluated 
with the current dataset. A survey of a larger area, and more fields at similar radius but 
different position angles would be necessary to ascertain whether it has been a global or a 
localized event. The fainter stars in the subgiant-branch region are well matched by the 13.5 
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Gyr, Z=0.001 isochrone, but no blue extended horizontal branch can be seen in the CMDs. 
The red clump — or red horizontal branch — presents, however, a short faint blue extension. 

Figure 3 displays the color function of the CMD integrated over the magnitude range 
21.5 < I < 17.5, i.e. including all stars from just below the red clump to the location of 
the oldest MS turnoffs. The star counts have been normalized to the number of stars in 
the last two bins, which mainly inchide the red clump and a portion of the RGB just below 
it. Since these structures contain stars of ages greater than ~ 1 Gyr, this is equivalent to a 
normalization by an indicator of the mean star formation rate. The gradient in the number 
of blue stars present in these CMDs, and their color distribution is very striking, notably 
more pronounced than would be expected from a visual inspection of the CMDs alone. The 
main population of age > 2.5 Gyr present in the GMD appears at {V — I) ~ 0.2, and 
the non-neghgible number of stars between {V — I) — — 0.2 indicate the presence of the 
population of age ~ 1.5 Gyr mentioned above (error-bars have been omitted for clarity, but 
they are smaller than the differences between the curves). Stars bluer than (V — I) ~ 0.5 are 
younger than ~ 8 Gyr, and it is in this interval of color that the largest differences between 
fields are apparent. Beyond R~470' the population of stars younger than 8 Gyr starts to 
abruptly decrease. The interval {V — I) = 0.6 — 0.7 is dominated by stars older than 8 Gyr, 
for which the relative numbers are very similar at all radii. Differences in the color function 
appear again in the interval {V — I) — 0.8 — 0.9, which mainly contains the subgiant branches 
of populations of all ages. 



5. Discussion 

We have presented a high-quality CMD, reaching the oldest main sequence turnoffs, 
of an LMC field situated at R=450~490' from the center of mass of the LMC, toward the 
North (i.e., in the direction of its major axis). Very few CMDs of the LMC field have been 
pubhshed beyond this radius. The most targeted, detailed work is, to our knowledge, that 
of Stryker (1984) who did photographic photometry of stars in an area of 118(4 x 118(4 
situated 9° northeast of the LMC bar. She just reached the oldest MS turnoffs. In spite of 
the very large photometric errors, she was able to show that a blue population brighter than 
the old turnoff was indeed present. In the very large area surveyed, she was also unable to 
find evidence of a blue horizontal branch. 

In this paper, we show the presence of a substantial intermediate- age stellar population 
at R = 8° (~ 7 kpc) from the center of the LMC. A strong gradient in this population, 
however, occurs in a range of just 0.5 kpc. The combination of an exponential surface 
brightness profile together with intermediate-age stars indicates the presence of a dominant 
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disk population even at these large galactocentric distances. It would be of great interest 
to investigate the stellar populations at even larger radial distances, to ascertain whether 
eventually a background old population appears and stays for a while (possibly indicating 
a true old stellar halo population that would need to be investigated further), or whether 
the actual outer edge of the LMC is smoothly reached with a continuing gradient of stellar 
populations. The lack of a well populated, blue extended horizontal branch, however, both 
in the field investigated here and in the larger, even more peripheral field studied by Stryker 
(1984) seems to indicate that the amount of field stellar population as old and metal poor 
as that of the Milky Way halo globular clusters and dSph galaxies is small in the LMC. 

Turning now to the origin of the intermediate-age population in the outer fields, it 
is important to ask whether it could have formed at its current position, or it must have 
migrated from a more internal birth site (Hodge, Hunter & Oey 2003). Since the rotation of 
the galaxy may make it difficult for stars to move outwards, this population of stars may be 
a sign that the HI gas extended as far as this outer field as recently as a few Gyr ago and has 
since (gradually, thus the observed gradient), receded to the more central regions. Regarding 
the youngest ~ 1.5 Gyr old population, which seems to be quite uniformly distributed over 
the field, an interaction between the LMC and the SMC, or the merger of the LMC with a 
smaller galaxy are attractive mechanisms for either in-situ formation, or accretion of stars 
in an apparently discrete age range. 

Our results, in any case, are a demonstration that the lack of MS stars and the presence 
of a RGB with characteristics very similar to that of an old population cannot be taken as a 
proof of the existence of an exclusively old stellar population. The presence of a true stellar 
halo in dirr galaxies remains therefore unproven, and it needs to be investigated with much 
deeper CMDs aided with spectroscopic data to characterize robustly the metaUicities and 
the kinematics of the individual stars. 

This research is part of a joint project between the Universidad de Chile and Yale 
University that is partially funded by the Fundacion Andes. C.G. acknowledges partial 
support from the lAC (Project 311902), the Spanish Ministry of Science and Technology 

(Plan Nacional de Investigacion Cientifica, DesarroUo e Investigacion Tecnologica, AYA2002- 
01939) and from the European Structural Funds. R.Z. acknowledges the support of NSF 
grant AST-0098428. 



-9- 



REFERENCES 

Aparicio, A., Tikhonov, N., & Karachentsev, I. 2000, AJ, 119, 177 
Aparicio, A. & Tikhonov, N. 2000, AJ, 119, 2183 
Bothun, G.D. & Thompson, LB. 1988, AJ, 96, 1988 

Gallart, C, Stetson, P.B., Pont, F., Zinn, R. & Hardy, E. 2004, in preparation 
Hardy, E. 1978, PASP, 90, 132 

Hidalgo, S., Marin- Franch, A. & Aparicio, A. 2003, AJ, 125, 1247 
Hodge, P.W., Hunter, D. & Oey S. 2003, PASP, 115, 273 

Irwin, M.J. 1991 in The Magellanic Clouds. lAU Symposium 148, 453 Ed. R. Haynes and 
D. Milne. Kluwer Academic Publishers, Dordrecht. 

Kim, S., Staveley-Smith, L., Dopita, M.A., Freeman, K.C., Sault, R.J., Kesteven, M.J., 
McConnell, D. 1998, ApJ, 503, 674 

Kunkel, W.E., Irwin M.J. & Demers, S. 1997, A&AS, 122, 463 

Landolt, A.U. 1992, AJ, 104, 340 

Letarte, B., Demers, S., BattineUi, P., Kunkel, W.E. 2002, AJ, 123, 832 
Minniti,D. & Zijlstra, A. 1996, ApJ, 467, L13 

Minniti,D. Borissova, J., Rejkuba, M., Alves, D.R., Cook, K.H. & Freeman, K.C. 2003, Sci, 
301, 1508 

Olszewski, E.W., Harris, H.C., Schommer, R.A. & Caterna, R.W. 1988, AJ, 95, 84 
Pietrinferni, A., Cassisi, S., Salaris, M. & Castelli, F. 2004,ApJ, in press 
Robin, A.C., Reyle, C, Derriere S. & Picaud S. 2003, A&A, 409, 523 
Sandage, A. 1971, ApJ, 166, 13 

Schommer, R., Smith, C, Olsen, K. & Walker, A. Eds. 2000: CTIO MOSAIC II Imager 
User Manual (http://www.ctio.noao.edu/mosaic/manual/index.html) 

Stetson, P.B. 1987, PASP, 99, 191 



stetson, P.B. 1994, PASP, 106, 250 
Stetson, P.B. 2000, PASP, 112, 925 
Stryker, L.L. 1984, ApJSS, 55, 127 

van der Marel, R. P., Alves, D. R., Hardy, E., SuntzefT, N. B. 2002, AJ, 124, 2639 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 11 - 




Fig. 1. — V (dots) surface brightness profile of the LMC as measured from our MOSAIC 
data. The field discussed in this paper is the outermost one. The surface brightness has been 
calculated across this field and given for seven sections at increasing distances from the center 
of the LMC. Two measurements of Hardy (1978, his innermost point and the outermost one 
of his "outer regions") have been added as triangles. The outermost radius measured by 
Bothun & Thompson (1988) is marked with an arrow, which also indicates an estimate 
of the V surface brightness measured by them at this radius, assuming {B — V) = 0.5. 
Exponential laws for the LMC disk as calculated from the whole dataset represented in the 
figure (Hardy 1978 plus our data, solid line; /ij(0)=21.16, a=98'.7) and as derived by Hardy 
(1978, long-dashed line; ;Uj(0)=21.18, a=86'2) and Bothun & Thompson (1988, short-dashed 
line; /ij(0)=21.16, a=100f6) are represented. 
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Fig. 2. — Observed [(V-I), I] CMDs in order of increasing distance from the LMC center. 
The number of stars in each CMD is the same and, since the density of LMC stars de- 
creases as a function of radius, 6R varies for each CMD. The range in R of each CMD is 
labeled. Isochrones by Pietrinferni et al. (2004) for Z=0.001 and ages 13.5 and 8 Gyr (red), 
Z=0.002 and age 4.5 Gyr (green) and Z=0.004 and ages 2.5 and 1.5 Gyr (blue) have been 
superimposed. 
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Fig. 3. — Color function of the CMD integrated over the magnitude range 21.5 < I < 17.5, 
for each of the CMDs shown in Figure 2. Sohd, long-dashed, short-long dashed, short- 
dashed, dot-short dashed and dotted lines are in order of increasing galactocentric distance, 
respectively (see text for details). 



